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SUMMARY 

This study develops an original and innovative matrix representation with respect to the information flow 
for networked multi-agent system. To begin with, the general concepts of the edge Laplacian of digraph are 
proposed with its algebraic properties. Benefit from this novel graph-theoretic tool, we can build a bridge 
between the consensus problem and the edge agreement problem; we also show that the edge Laplacian 
sheds a new light on solving the leaderless consensus problem. Based on the edge agreement framework, the 
technical challenges caused by unknown but bounded disturbances and inherently nonlinear dynamics can 
be well handled. In particular, we design an integrated procedure for a new robust consensus protocol that is 
based on a blend of algebraic graph theory and the newly developed cyclic-small-gain theorem. Besides, to 
highlight the intricate relationship between the original graph and cyclic-small-gain theorem, the concept of 
edge-interconnection graph is introduced for the first time. Finally, simulation results are provided to verify 
the theoretical analysis. Copyright © 2010 John Wiley & Sons, Ltd. 
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1. INTRODUCTION 

Recent decades have witnessed tremendous interest in investigating the distributed coordination 
of multi-agent systems due to broad applications in several areas including formation flight [2], 
coordinated robotics [9], sensor fusion [18] and distributed computations [17]. Recently, the 
consensus, as a critical research problem, has received increasing amounts of attention. To address 
such issue, the study of individual dynamics, communication topologies and distributed controls 
plays an important role. 

Graph theory contributes significantly to the analysis and synthesis of networked multi-agent 
systems because it provides general abstractions for how information is shared between agents in 
a network. Particularly, the graph Laplacian plays an important role in the convergence analysis 
and has been explored in several contexts [19][21]. Although the graph Laplacian is a convenient 
method to describe the geometric interconnection of networked agents, another attractive notion- 
the edge agreement, which has not been explored extensively, deserves additional attention because 
the edges are adopted as natural interpretations of the information flow. Pioneering studies on edge 
agreement protocol provide new insights into how certain subgraphs, such as spanning trees and 
cycles, affect the convergence properties and set up a novel systematic framework for analyzing 
multi-agent systems from the edge perspective [31][32][34]. In these literatures, an important matrix 
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representation is introduced, which is referred to as edge Laplacian. The edge agreement in [31] 
provides a theoretical analysis of the system’s performance using both H 2 and H ^ norms, and 
these results are applied in relative sensing networks, referring to [32]. Furthermore, based on 
the properties of the edge Laplacian, [34] examines how cycles impact the H 2 performance and 
proposes an optimal strategy for designing consensus networks. It is worth noting that, in above- 
mentioned studies, the edge Laplacian representation is valid for undirected graphs, rather than 
the more general directed case. Considering the numerous applications of digraphs in multi-agent 
coordination, extending such a graph tool to the directed case is of great interest. 

The information exchange between agents cannot be accurate in practical applications because 
it is inevitable that types of external noises exist in sensing and transmission; and the existence of 
perturbations always leads to instability and performance deterioration. Therefore, it is significant 
to investigate their effects on the behavior of multi-agent systems and design a robust consensus 
protocol to improve disturbance rejection properties. In [1], the concept of e-consensus for an 
undirected graph is introduced in which the state errors are required to converge in a small region 
under unknown but bounded external disturbances. In the presence of external disturbances, [29] 
guarantees a finite L 2 -gain performance under strongly connected graphs. The global consensus 
with guaranteed lL Xj performance can be achieved under a strongly connected graph in [8], where 
a Lipschitz continuous multi-agent system subject to external disturbances is considered. Motivated 
by above observations, we consider a more general directed communication graph containing a 
spanning tree and has unknown but bounded disturbances in the neighbors state feedback. 

Recently, researchers have increasingly interested in consensus for multi-agent systems with 
nonlinear dynamics because most of the physical systems are inherently nonlinear in nature. In [4], 
an adaptive pinning control method is introduced to study the synchronization of uncertain nonlinear 
networked systems. The second-order consensus of multi-agent systems with heterogeneous 
nonlinear dynamics and time-varying delays is investigated in [36]. In [27], a distributed control 
law for the formation tracking of a multi-agent system that is governed by locally Lipschitz- 
continuous dynamics under a directed topology is developed. It should be noted that the studies 
mentioned above are all in the leader-follower setting; however, we will investigate a substantially 
challenging problem in this paper, i.e., the leaderless consensus problem for multi-agent systems 
with inherently nonlinear dynamics. To the best of our knowledge, most of existing studies focus 
on solving the nonlinear leaderless consensus under undirected graphs [7][14][20], while a few 
considered the directed topology. Only a small number of studies use the concepts of generalized 
algebraic connectivity [30], nonsmooth analysis based on the study of Moreau [10] and the limit- 
set-based approach [22] to solve the state agreement problem with nonlinear multi-agent systems 
under digraphs. However, due to the fact of these methods’ extremely complicated characters, it is 
difficult to promote them in the practice. The edge agreement model allows the development of a 
universal solution framework to address these problems. 

To address the technical challenges caused by the external disturbances and the inherently 
nonlinear dynamics, the concept of input-to-state stability (ISS), which reveals how external inputs 
affect the internal stability of nonlinear systems (see [23] for a tutorial), is used in this paper. 
Actually, the concept of ISS is of significance in analyzing large-scale systems, and considerable 
efforts have been devoted to the interconnected ISS nonlinear systems; particularly, the general 
cyclic-small-gain theorem of ISS systems developed in [1 1 ]. Most recently, the nonlinear small-gain 
design methods, especially the cyclic-small-gain approach, are utilized to design new distributed 
control strategies for flocking and containment control in [28] as well as to deal with formation 
control of nonholonomic mobile robots in [12]. In [13], the authors present a cyclic-small-gain 
approach to distributed output-feedback control of nonlinear multi-agent systems. 

This paper initially extends the concept of the edge Laplacian to digraphs and explores its 
algebraic properties. To proceed with a seamless integration of graph theory and ISS designs, we 
present a new robust consensus protocol induced from the edge agreement model to solve the 
leaderless consensus problem with unknown but bounded disturbances and inherently nonlinear 
dynamics. To better comprehend the edge agreement mechanism, both strongly connected and 
quasi-strongly connected situations are considered in this paper. The contributions of this paper 
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depend on three aspects. First, the edge Laplacian of digraph is proposed as well as the edge 
adjacency matrix. Comparing to our previous works [35], much more details about the algebraic 
properties of the edge Laplacian are explored. As a matter of fact, the invertibility of the incidence 
matrix and the spectra properties of the edge Laplacian play a central role in the subsequent analysis. 
Second, we provide a general framework for analyzing the leaderless consensus problem for digraph 
based on the edge agreement mechanism. Under such framework, the technical challenges caused 
by the external disturbances and the inherently nonlinear dynamics can be effectively addressed. 
In particular, we design an integrated procedure for a new robust consensus protocol which is 
based on a blend of algebraic graph theory and ISS method. Third, following our setup, the edge- 
interconnection graph is proposed. One of our primarily goal in this paper is to explicitly highlight 
the insights that the edge-interconnection structure offers in the analysis and synthesis of multi¬ 
agent networks. In this direction, we note a reduced order modeling for the edge agreement in terms 
of the spanning tree, and based on this observation, a two-subsystem interconnection structure is 
given. 

The organization of this paper is as follows: In Section 2, a brief overview of the basic concepts 
and results in graph theory are presented as well as the ISS cyclic-small-gain theorem. The edge 
Laplacian of digraph and other related concepts are proposed in Section 3. The main content on the 
robust consensus protocol with inherently nonlinear dynamics is elaborated in Section 4. Numerical 
simulation results are provided in Section 5. The last section presents the conclusions and proposes 
a number of future research directions. 


2. BASIC CONCEPTS AND PRELIMINARY RESULTS 

In this section, we present a number of basic concepts in graph theory and the ISS cyclic-small-gain 
theorem. 


2.1. Graph and Matrix 


In this paper, we use |-| and ||-|| to denote the Euclidean norm and 2-norm for vectors and matrices 
respectively. The notation is used to denote the supremum norm for a function. Let 1Z ( A ) and 
M (A) denote the range space and null space of matrix A. 


Let Q = (V. £, Ag) be a digraph of order N with a finite nonempty set of nodes 

V = {1,2, • • ■ , N}, a set of directed edges £ C V x V with size L and a adjacency matrix 
Ag = [ay] e R NxN , where ay = 1 if and only if (i,j) € e else ay = 0. The degree matrix 


N 

_ ^7 = 


•j =1 ai J > 


. i = 1,2, • • • ,N, being the out-degree 


Ag = [Ay] is a diagonal matrix with [A* 
of node i, and the graph Laplacian of the digraph Q is defined by Lg = Ag — Ag. For 
a given parent node i, its incident edge is called the parent edge denoted by V\ its 
emergent edge is called the child edge denoted by C®; and the edges derived from the same 
parent i (i.e., the collection of C) are called the sibling edges denoted by £' c . We call 
an edge the neighbor of e& if they share a node, and the neighbor set of ty is denoted 
by N f:k = {ei £ £ : if cy and ci share a node}. The outgoing neighbors of ty refer to Nf = 
{ei e N ek : if efc = V 1 and ei = C®; or e*, ej £ £g, where i is the conjunct node} (see, e.g.. Figure 
1). 




Figure 1. Two typical forms of the outgoing neighbors, N^ k 
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A directed path in digraph Q is a sequence of directed edges. A directed tree is a digraph in which, 
for the root i and any other node j, there exists exactly one directed path from i to j. A spanning 
tree of a digraph is a directed tree formed by graph edges that connect all the nodes of the graph [5]. 
Graph Q is called strongly connected if and only if any two distinct nodes can be connected via a 
directed path and quasi-strongly connected if and only if it has a directed spanning tree [25]. 

The incidence matrix E (Q) for a digraph is a {0, ±l}-matrix with rows and columns indexed by 
the vertexes and edges of G, respectively, such that 


[E{Q)] ik = { 


+1 

-1 

0 


if i is the initial node of edge 
if i is the terminal node of edge 
otherwise 


which implies each column of E contains exactly two nonzero entries “+1” and “-1”. In addition, 
for a quasi-strongly connected digraph, the rank of the incidence matrix is rank(E ((?)) = N — 1 
from [25]. Figure 2 depicts an example with its incidence matrix. 



2.2. ISS and Small-gain Theorem 

A function a : M+ —> K + is said to be positive definite if it is continuous, a(0) = 0 and a(s) > 0 for 
s > 0. A function 7 : —> R + is of class /C if it is continuous, strictly increasing and 7 (0) = 0; it 

is of class /Coo if, in addition, it is unbounded. A function fi : R + x K + -» R + is of class ICC if, for 
each fixed t, the function 8 (•, t) is of class K. and for each fixed s, the function 8 (s, •) is decreasing 
and tends to zero at infinity. Id represents identify function, and symbol o denotes the composition 
between functions. 

Consider the following nonlinear system with x € R" as the state and w e K m as the external 
input: 

x = ((x,w) ( 1 ) 

where ( : M" x K m —y M n is a locally Lipschitz vector field. 

Definition 1 ([ 23 ]) 

The system (1) is said to be input-to-state stable (ISS) with w as an input if there exist 8 G /C£ and 
7 e /C such that for each initial condition (0 j and each measurable essentially bounded input «;(■) 
defined on [ 0 , 00 ), the solution x(-) exists on [ 0 , 00 ) and satisfies 

k(f)l < /3(k(0)|,7(lk||oo)), V/ > 0. 

It is known that if the system x = f(x,w) in (1) is ISS with w as the input, then the unforced 
system x = f(x, 0 ) is globally asymptotically stable at x = 0 . 

Lemma 1 ([ 23 ]) 

The system (1) is ISS if and only if it has an ISS-Lyapunov function. 

Consider the following interconnected system composed of N interacting subsystems: 

Xi = (i(x,Wi), i = 1,... ,N ( 2 ) 
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and C i ■ R n+mi 




where Xi € R ni , w, G K’ 

continuous such that x = [x[ 
condition. The external input w = 
function from R + to with m = YliLi 


with 


v N 

n = i m 


is locally Lipschitz 


is the unique solution of system (2) for a given initial 


function from ./.--subsystem to //-subsystem [6]. 


, w/^] T is a measurable and locally essentially bounded 
rrii. Furthermore, we use 7 * G 1C to represent the gain 


Lemma 2 (Cyclic-small-gain Theorem, [11]) 

Consider the continuous-time dynamical network (2). Suppose that for i = 1,... ,N, the x,- 
subsystem admits an ISS-Lyapunov function V, : IR"* —y R satisfying 


• there exist a i; a* € /Coo such that 


—i (I R- 1) — (Xi ) S ctj(|x^|), 


• there exist 7 G 1C U {0} (j ^ i), 7 ™/ G 1C U {0} and a positive definite on such that 

Vi(xi) > max{'yZi(V j (x j )),'Y™ i (\w i \)} 

=^VVi(xi)(i(x,Wi) < Vx, Vwj. 


Then, the system (2) is ISS if for each r = 2,... TV, 

X^o Xji t 1 

7x 4l o 7x i2 0 • • • 0 7x ir < Id 


(3) 


for all 1 < ij < N, ij ^ ij< if j ^ j'. 

3. THE EDGE LAPLACIAN OF DIGRAPH 

The edge Laplacian is a promising graph-theoretic tool, however it still remains to an undirected 
notion and is thus inadequate to handle our problem. Undoubtedly, extending the concept of the 
edge Laplacian to the digraph and exploring its algebraic properties will contribute significantly 
to the investigation of multi-agent systems. In this section, we first give the definition of the in¬ 
incidence matrix and out-incidence matrix. 

Remark 1 

As known that the “out-degree” relates how each node in the network impacts on other nodes, but 
the “in-degree” directly captures how the dynamics of an agent is influenced by others [15]. In fact, 
the following investigations, including the definitions and properties, all adopt the “out-degree” 
description. However, for the “in-degree” case, analogous methods can also be applied. 

Definition 2 ( In-incidence and Out-incidence Matrix) 

The N x L in-incidence matrix Eq (G) for a digraph Q is a {0,-1} matrix with rows and columns 
indexed by nodes and edges of G, respectively, such that 

J — 1 if i is the terminal node of edge efc 
[E0(5)1 “ := \o otherwise 

and the out-incidence matrix is a {0, +1} matrix, defined as 

j+1 if z is the initial node of edge efc 
[ ® <5)1 “ := \o otherwise. 

Comparing with the definition of the incidence matrix, we can write E (Q) in the following 
manner 


E(G) = Eq(G) + E®(Q). 


(4) 
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In the following discussions, we use E , Eq and E® instead of E(Q), E@(G) and E®(Q). 

Due to the fact that each row of the out-incidence matrix actually can be viewed as a 
decomposition of the out-degree from a node to each specific edge, we can derive a novel 
factorization of the graph Laplacian. 

Lemma 3 

Considering a digraph Q with the incidence matrix E and the out-incidence matrix E®, then the 
graph Laplacian of Q have the following expression: 


Lg — E^E T . 


(5) 


Proof 

According to the preceding definition of E<%, we obtain that 




T 


[Aii] if i = j 

0 otherwise 


which implies E®E® = Ag. Similarly, we have 


E, 




if j£Ni 

otherwise 


where N, denotes the neighbor set of node i, and we can collect the terms as E®Eq = —Ag. 
According to the definition of the graph Laplacian and equation (4), we have 

E®E T = EqEq + EqEq = Ag — Ag = Lg. 


The proof is concluded. □ 

Next, we give the definition of the edge variant of the graph Laplacian. 

Definition 3 (Edge Laplacian of Digraph) 

The edge Laplacian of digraph is defined as 


L e E t E® 


with L x L elements. 

Definition 4 (Edge Adjacency Matrix ) 


+1 



-1 

0 


e k ,ei G £' c 

e k = V 1 and e; = C l 

otherwise. 


For instance, the edge Laplacian matrix and the edge adjacency matrix of the simple digraph 
shown in Figure 2 are, respectively 


f 1 

-1 

1 



(° 

-1 

1 

_1 \ 

-1 

1 

-1 

1 

and A e = 

-1 

0 

-1 

1 

1 

0 

1 

0 

1 

0 

0 

0 

V 0 

1 

0 

l) 


V 0 

1 

0 

0 / 


Lemma 4 

The edge Laplacian can be constructed from the edge adjacency matrix 


L e — I + A e . 


( 6 ) 
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Proof 

According to the definition of 77® and A e , the result comes after the proof of Lemma 3. □ 

To provide a deeper insight into what the edge Laplacian L e offers in the analysis and synthesis 
of multi-agent systems, we propose the following lemma. 

Lemma 5 

For any digraph Q, the graph Laplacian Lg and the edge Laplacian L e have the same nonzero 
eigenvalues. In addition, the edge Laplacian L e contains exactly TV — 1 nonzero eigenvalues and all 
in the open right-half plane, if Q is quasi-strongly connected. 

Proof 

Suppose that A ^ 0 is an eigenvalue of Lg, which is associated with a nonzero eigenvector p. 
Therefore, we have 


Lgp = E®E T p = A p (7) 

which implies E T p = p ^ 0. By left-multiplying both sides of (7) by E T , one can obtain 

L e p = E T E®E T p = A p 

which shows that L e contains the nonzero eigenvalues that Lg has. 

By using similar approaches, we can proof that Lg also has all the nonzero eigenvalues of L e . It 
turns out that the nonzero eigenvalues of Lg and L e are identical. 

By Lemma 3.3 in [21], for a quasi-strongly connected digraph Q of order TV, Lg has TV — 1 
nonzero eigenvalues and all in the open right-half plane, therefore L e contains exactly TV — 1 
nonzero eigenvalues as well. Then we come to the conclusion. □ 

Lemma 6 

Consider a quasi-strongly connected digraph Q of order TV, the edge Laplacian L e has L — TV + 1 
zero eigenvalues and zero is a simple root of the minimal polynomial of L e . 

Proof 

Consider the quasi-strongly connected digraph Q , L e has exactly TV — 1 nonzero eigenvalues; 
therefore 

rank(L e ) > TV — 1 (8) 

and the algebraic multiplicity of the zero eigenvalue of L e is L — TV + 1. Besides, recall the fact that 

L e = E t E ® and rank(E) = TV — 1, so we have 


rank(L e ) < rank(E T ) = TV — 1. (9) 

Then by combining (8) and (9), one can obtain rank(L e ) = TV — 1, i.e., the dimension of the 
null space of L e is dimAf(L e ) = L — TV + 1. In other words, the geometric multiplicity of zero 
eigenvalue is L — TV + 1. Clearly, the geometric multiplicity and the algebraic multiplicity of zero 
eigenvalue are equal, which implies that the corresponding Jordan block for each zero eigenvalue is 
size one from [16]. That is, zero is a simple root of the minimal polynomial of L e . □ 


From [33], for the undirected graph, the dynamics of the edge agreement model can be captured 
by the reduced order system based on the spanning tree. Next, we will further discuss the similar 
results under digraph. 

Clearly, if the digraph Q is quasi-strongly connected, it can be rewritten as a union form: 
Q = QpU Qc, where Qj- is a given spanning tree and Qc is the cospanning tree respectively. 
Correspondingly, the incidence matrix can be rewritten as 

E = [ Ej- Eg] 
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through some permutations, and Ex, Ec are incidence matrices with respect to Qj- and Qc. Similarly, 
the out-incidence matrix can be rewritten as 

-E® = [E® 7 - E®c] . 

It should be mentioned that Qc can be reconstructed from Qp, which indicates that Ej- has full 
column rank [33]. 

According to the partition, one can represent the edge Laplacian in terms of the block form of the 
incidence matrix as 


E£E® r E^E®c 


L e l L e 2 

E^ E® 7- E^ E®c. 


Le 3 L e 4 


Also it is useful to express the edge adjacency matrix A e as the block representation 


a _ A e i A e 2 
e _ \_A e3 A e4 ' 

Following from ( 6 ), we have 

L e i = E + A e \, L e 2 = A e 2 , L e 3 = A e 3 , and L e4 = I 2 + A e4 . 

Lemma 7 

Considering a quasi-strongly connected digraph Q = Qj- U Qc, the pseudoinverse of the incidence 
matrix if exists, and there exists a matrix R such that 

E^ = R^E\- (10) 

with l A = if [Rif] 1 and lf r = [R^If] 1 lf r , where if is the right-inverse of R and lf r is 
the left-inverse of E 7 -. 


Proof 


Since Ej- has full column rank, so its left-inverse exists and can be directly obtained by 


Er = [Ej-Er]- 1 Ej-. Because the columns of Ec are linearly dependent on the columns of E 7 -, 


we have 


E t T = E c (11) 

and then the matrix T can be obtained by 

T = E|E C . (12) 


The matrix R 
cut space of Q 

is now defined as 
[5]. Besides, the 

R = [I T \. In fact, the rows of the matrix R form a basis for the 
incidence matrix of Q can be written as E = JfR. Clearly, Ej- is 

of full column rank and R is of full row rank; therefore, the pseudoinverse of E can be calculated 


by Et = R^E\- from [3]. Then we reach the conclusion. 


□ 


4. ROBUST CONSENSUS OF NONLINEAR MULTI-AGENT SYSTEMS VIA ISS DESIGN 

In this section, a new consensus protocol is presented through a seamless integration of graph 
theory and ISS design. The newly developed cyclic-small-gain theorem is employed to address the 
challenges caused by unknown but bounded disturbances and the inherently nonlinear dynamics. 
Contrary to the well-studied graph Laplacian dynamics, we will analyze and synthesize multi¬ 
agent systems with the edge perspective by using edge agreement framework. To facilitate a better 
understanding, both strongly connected and quasi-strongly connected situations are considered. 
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To begin our analysis, the dynamics of the i-th agent is defined as 


Xi ft) = f (t, Xi) + Wi (t) + fii(t),i = 1,2, • • • ,7V (13) 

where Xi ( t ) e t n refers to the state vector of the *-th node; / (t,Xi) : 1 x 1" -) R n denotes a 
Lipschitz continuous function; (t) e R" describes unknown but bounded disturbances with the 
upper bound £ > 0, i.e., | Wi (f)| < £; m (t) £ R" represents the control input. 

In the presence of the disturbances, we should not expect that agents can accurately reach 
consensus. Therefore, we introduce the robust consensus to describe the influence of the 
disturbances on the behavior of the system. 

Definition 5 (Robust Consensus ) 

In the presence of disturbances w(t) with the upper bound £, we design a distributed control law 
Hi ( t ) for i = 1,2, ■ • • ,7V such that the agent state x, governed by (13) can reach robust consensus 
in the nonlinear-gain sense that 

kift) ~Xj(t)\ <^(|K*)|| 00 ), t CK, as Moo. 

Before proceeding, we make the following assumption. 

Assumption 1 

For the nonlinear function / (i, Xi) in (13), there exists a nonnegative constant rj such that 
1/ ft, xi) - f(t,x 2 )| <r]\xi- x 2 \ ,Vxi,x 2 e R n ;t > 0. 

To reach consensus, we employ the following distributed consensus protocol: 

Mi ft) = - {Xi ft) - Xj ( t )) +Ui ft),* = 1,2,--- ,7V (14) 

jCNi 

where Ui (t) is the auxiliary control input yet to be designed. 

4.1. Edge Agreement 

Considering an edge e*,, we define the edge state as x f . k ( t ), which represents the difference between 
two agents associated with e*,. We then have 

£e fe ft) = Z®(e fc) ft) - a; 0( e fc) ft) (15) 

where Cg) (efc) and 0 (e*,) denote the initial node and the terminal node of e*,, respectively. 
Following this, we obtain 

x e {t) = E T x{t) (16) 

where x e (t) is the collection of x ek ( t ). Consider the well-known graph Laplacian dynamics (linear 
and noise-free) in [19] as 


x ( t ) = —Lx (t). (17) 

Differentiating (16) and substituting in (17), leads to 

x e ft) = -L e x e (t) (18) 

which is referred to the edge agreement protocol. In comparison to the consensus problem, the edge 
agreement, rather than requiring the convergence to the agreement subspace [19], expects the edge 
dynamics (18) to converge to the origin. Essentially, the evolution of an edge state depends on its 
current state and the states of its adjacency edges. In addition, the edge agreement of x e implies 
consensus if the digraph Q has a spanning tree [31]. 
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Remark 2 

Recall that the objective leaderless consensus is defined as lim \xi (t) — x :j (i)| = 0. Obviously, 

t—> OO 

directly modeling such problems is difficult since the coupling of different agents Xi (t) and x :j (t) 
is involved. However, the edge agreement framework provides a possible method of studying the 
leaderless consensus problem from the edge perspective. By using (15) and (16), we can turn the 
consensus problem into an edge agreement problem, where the asymptotic stability of x e implies 
the consensus. In fact, based on such framework, the leaderless consensus problem can be extremely 
simplified. Moreover, we also suggest that the edge agreement does not only provide a broader scope 
for addressing the leaderless consensus problem but also has a significant potential to address the 
leader-follower case. 

Given the protocol (14), the graph Laplacian dynamics is obtained as 

X (t) = T ( t, x) — LgX ( t) + w (t) + u ( t ) (19) 

where T {t,x), w (t) and ^ (t) are the column stack vectors of f(t,Xi), Wi(t ) and for 

*= 1,2,--- ,N. 

By differentiating (16) and substituting into (19), we have the following edge Laplacian dynamics: 

x e ( t) = E T T (t, x) — L e x e (t) + E t w (t) + E t u ( t) (20) 

where u (t) needs to be further determined. For convenience, we define 

u e ( t) = E t u ( t) . (21) 


From (6), we have 


Le k — 1 + ^ [A e ]kl- 

e,€ N ® k 

Finally, the sate of the fc-th edge evolves according to the following system: 


Xe k (t) = f (t, X® (ek) ) - / (f, X Q{ek) ) - j X ek (t) + ^ [A e ] kl X e , (t) + [E T ] k w(t)+ U £k (t) . 

e l6 JV® J 

( 22 ) 



Figure 3. The corresponding edge-interconnection digraph of Figure 2. 


Remark 3 

We use fx (t) as the only implementable control input for each agent, while use u e in the following 
analysis as a matter of convenience. However, we cannot presume the existence of u (t) because 
(21) may not have a solution. Actually, whether or not the equation has a solution depends on the 
topological structure of the graph, and details will be discussed in the next section. 

By translating into the edge Laplacian dynamics, the novel system establishes a new control 
interconnection relation. To describe this relationship, we provide the following definition. 
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Definition 6 ( edge-interconnection digraph ) 

Considering the x ek -subsystems as nodes and the control interconnections as directed edges, the 
interconnected system composed of the x eic -subsystems can be modeled as a digraph Q called as 
edge-interconnection digraph. 

In fact, Q can be easily constructed from Q, and equation (22) illustrates the connection of the 
edge-interconnection digraph. In particular, there are two steps to transform Q into Q. First, all 
the edges e* in Q are modeled as the vertices and denoted by tt*.. Second, base on (22) and the 
definition of A e , for any specified edge pair e^, e/, if eh, ei € ££, then and tt/ are connected by a 
bidirectional edge; and if e& = V and e; = C\ there will be a directed edge incident into tt/, from tt/. 
Figure 3 depicts the corresponding edge-interconnection digraph of the example shown in Figure 2. 
To reveal the intricate relation between the original digraph and the edge-interconnection digraph, 
we need the following Lemma. 

Lemma 8 

For a rooted tree, the corresponding edge-interconnection digraph Q is composed of several strongly 
connected subgraphs Q, , i = 1,2, • • •, which are coupled through simple cascaded connections and 
parallel connections. 

Proof 

Obviously, the sibling-edges E' c with the same parent i form a strongly connected components in the 
edge-interconnection digraph, since they affect each other mutually. Besides, for the e& = V 1 , and 
d = C\ there will be a directed edge incident into tt/, : from tt/ which forms the cascaded connection 
in the same branch. When taking the strongly connected components as nodes, Q is acyclic and 
consists of several cascaded connections and parallel connections. For instance, consider a rooted 
tree in Figure 4(a), the corresponding edge-interconnection digraph is illustrated in Figure 4(b). □ 



Figure 4. The corresponding edge-interconnection digraph with 3 strongly connected components , Q 2 

and £ 3 . 


Lemma 9 

The edge-interconnection graph Q is strongly connected if and only if Q is strongly connected 
digraph. 


Proof 

Clearly, for the strongly connected graph, there will be a directed path connecting any pair of edges 
(e/c, ej) in each direction. In that way, any two distinct nodes ( 7 ^, 7 q) of Q can be connected via a 
directed path; therefore, Q is strongly connected. On the other hand, while Q is strongly connected, 


it implies that, for any pair of nodes of Q, there always exists a directed path connecting them, i.e., 
Q is also strongly connected. □ 
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4.2. Main Results 

In this section, both strongly connected digraphs and quasi-strongly connected digraphs are 
considered. The cyclic-small-gain theorem is then employed to guarantee the robust consensus 
of the closed-loop multi-agent systems. Note that ISS cyclic-small-gain theorem can be directly 
applied if the underlying digraph is strongly connected from [11]. However, for the quasi-strongly 
case, we will translate it into a two-subsystem interconnection structure. 

4.2.1. Strongly Connected Digraph A digraph Q of a six-agent system is shown as an example 
in Figure 5(a), and the corresponding edge-interconnection digraph is shown in Figure 5(b). From 
Lemma 9, we know that the edge-interconnection digraph is strongly connected if Q is strongly 
connected. 



(a) The strongly connected digraph. 


Figure 5. The original strongly connected digraph 



(b) The edge-interconnection digraph, 
the corresponding edge-interconnection digraph. 


To begin our analysis, we define the following ISS-Lyapunov function candidate 

1 _ T _ 

= ) k = 1,2, ■ ■ ■ , L 

and we denote 0 as the set of all simple loops (more details please refers to [1 1 ]) of Q, and A a ( 7 -'*) 

e fc 

as the product of the gain assigned to the edges of a simple loop o £ 0, where 7 ~ e ' £ with 
fc = 1,2, ■ ■ ■ ,L and ej £ IV®. 

The main result for the strongly connected digraph is given as follows. 

Theorem 1 

Assuming that the digraph is strongly connected, consider the subsystem (22) with x ek as the state, 
and x ei (ei £ iV®) as the external inputs. For any specified constant a ek and 77 (ej £ IV®), 
we can design 


v\x ek \+ UeUpZ] ^\) + ^ 

ei£N®, 


+ 1 - 




(23) 


-1 


with p-‘ k = a 

x h — 

function V ek satisfying 


^ 7 j°‘ ^ oa(s) such that the subsystem (22) is ISS with ISS-Lyapunov 

{•AlM 


V e , > max 

e,eJV e ®. 
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If the required cyclic-small-gain condition (3) is satisfied, the composed system is ISS. Then, the 
objective robust consensus can be achieved by using the distributed consensus protocol ( 14) with 

u =[E t ] V (24) 


Proof 


Since V ek 


> max 

eieN®, 


R «.)}■ 


we have 


l*e.l < Pill ' ei e N Z- 

Note that | [E T ] | = y/2 and by taking the derivative of V ek , we have 

VI4 fe X efc — X ek Xe k 


= Xe h [Ue„ + (/ x ®(e k )) ~ / (*> ®0(e fc ))) + [ E T } k W - ( X ek + E [R]fe^e, ] ] 

ei£N® 


< x e k (lie* -®ej + r? 


RJ \ x e k \ + 




E tM. 




i£N® 


<i T e k {u ek -x ek )+-q\xl k \\x ek \ + V2t\x T ek \ + \xl k \ E IRl/cilPsI" (I^J) • 

ei£N® 


Using (23), we have 


VlRx efc S ^ X e k x e k ~ a e k ^e k 

which implies that the x^-subsystem is ISS. 

Since the induced edge-interconnection digraph is strongly connected, the ISS cyclic-small-gain 
theorem can be directly implemented. If the following cyclic-small-gain condition is satisfied 

MiZ l k ) < Id 


then the composed system (20) is ISS. 

It should be mentioned that as u e is designed, the whole system x e (t) is unforced. Based on the 
ISS property, we have 

\x e (t )| < /3(|x°| ,t) + 7 (MR),* > 0 

where x° is the initial state of x e (t) and j3 £ ICC, 7 £ 1C. Obviously, as t —> 00 , we have 
/3 (|x°| , t) =0. So that 


lim \x e {t)\ <7(IMIoo)> 

t—foo 

which implies the robust consensus. Since Q has a spanning tree, the pseudoinverse [E T ]^ 
exists. Then we can obtain the implemented consensus control input by using (24). The proof is 
concluded. □ 

4.2.2. Quasi-strongly Connected Digraph In this section, we consider a quasi-strongly connected 
digraph Q = Gr U Qc- An example is given as Figure 6 (a). The edges of Gt are marked as red. 
Accordingly, the edge-interconnection digraph is shown in Figure 6 (b), which consists of two parts: 
Q = Q r u Qc- Correspondingly, the edge Laplacian dynamics system x e (t) can be modeled as the 
interconnection of the Hj-- subsystem and the //c-subsystcm based on Gr and Qc- 
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(a) The quasi-strongly connected (b) The edge-interconnection 

digraph. digraph. 


Figure 6. The quasi-strongly connected digraph and the corresponding edge-interconnection digraph. 


As previously mentioned, the incidence matrix can be rewritten as E = [ Ej- Eq\ , and the edge 
Laplacian can be represented as the block form L e = r el e2 in line with the permutation. 

^e3 -^e4 

Therefore, the edge Laplacian dynamics x e (t) described by (20) can be translated into the following 
form: 


H T : x r (t) = T r ( t , x) - L ei x r (f) - L e2 x c ( t ) + E^w ( t ) + u r ( t ) 

with Tr (f, x) = Ej-E ( t , x), ur (f) = E^-u ( t ). 

He : xc ( t) = JF C (t, x) - L ei x c (t) - L ea x r (t) + E^w{t) + uc (t) 

with Tc (t,x) = E?j-E(t,x), uc (t) = E^u ( t ). Besides, E^w and E£w indicate unknown but 
bounded disturbances on Qr and Gc, respectively. The interacting system is shown in Figure 7. 






y*c 

• Xr 


fj 

/ x r 


H n 


e 


+ u c 


Figure 7. The interconnection of subsystems Ht and H c . 


Obviously, for a specific spanning tree, it contains TV — 1 edges. For subsystem xr k , with k = 
1 ,2 ,-•• ,TV - 1, ei,e k £ Gr, we have 


x r k = F Tk 


x Tk 


+ ^2 l A ei]kEE] ~ Le 2k xc+[E^\ k w + ur k . 


eieK 


We choose the following ISS-Lyapunov candidates: 

= ij x r k x r k , k = 1,2, - ■ ■ , TV — 1 
1 T 

v t = 7; x r x r 

1 T 

Vc = 2 X c x C - 


(25) 


(26) 

(27) 

(28) 
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Denote © 7 - as the set of all the simple loops of Gy, and denote Ay o ( 7 -) as the product of the 
gain assigned to the edges of a simple loop in ©7-. Let Ai denotes the smallest nonzero eigenvalue 
of matrix TT t , where T is defined in (12). 

We then present the main result for the quasi-strongly connected digraph as follows: 

Theorem 2 

Assuming that the digraph is quasi-strongly connected, consider the subsystem (25) with 
xy k (e/t € Gr) as the internal state. Let xy (e; £ Nf k fl G r ) and xq (e c £ Gc) be the external inputs. 

For any specified constant cry k > 0 and 7 - 71 , 7 ?° £ k„o, we can design 

' k ' k 


u T k = 


%T k 

\x%\ 


V\XT k \+ Peilfc/I^r? (\xT k \) + \Le 2k \pZ k (\XT k \)+V2£ 


ei£N 


® 

e k 



a r k 

2 


) 


XT k 


(29) 


with = a 1 o 0 a ( s ) anc * p7c = — 1 0 (^r fc ) 0 01 ( s )» suc ^ subsystem 

(25) is ISS with an ISS-Lyapunov function Vy k satisfying 


Vy k > max 

ei£N® 


7 11 ( Vr ,),< 


Ofc)} 


v v r k xr k < -VT k Vr k - 


If the cyclic-small-gain condition (3) as well as 


lit < 


XT * Ai (N - 1) 


(30) 


are satisfied, the whole system is ISS. Then the objective robust consensus can be achieved by using 
the protocol (14) with 

u =[£?] t « r . (31) 


Proof 

The main proof procedure contains three steps. Firstly, for the xy k -subsystem, the ISS properties 
can be guaranteed by taking (29) as the control law. Second, the ISS properties of the upside 
subsystem H-y are then proven by utilizing the ISS cyclic-small-gain theorem. Finally, we prove 
that the downside subsystem If: is ISS, and the objective robust consensus can be achieved while 
the small gain condition 7 ? c o -/If < Id is satisfied. 

Step 1 : 

Using the Lyapunov candidate defined in (26) and considering Vy k > 

max 1 7^' (V77), H c Tk (U c )| with e h e k £ Gr, e c £ Gc , we have 

e i€-N® k ^ ^ 


\xc\< P Z k (\xr k \) 
\xr,\ < p7y (\xr k \) ■ 
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Taking the derivative of Vp k , we have 


VVr fc ir fc = xp k x Tk 


~T 
= x T k 


u T k + f (t, X®(e fc) ) - / {t, x o(e k )) ~ \ Xp k + X l A ^kl x Ti ~ Le 2 k X C + [E%] } 


VJ 


eiGN ? 


< ^ ( u r fc - *rj + V \xp k | I xp k \~xp k X - x^- k Le 2 k x c + x\ [E %], 


e,£JV| 


< x r k ( u r k ~xr k ) + v\ x r k \ \ x T k \ + |*rj XI PeJJp^ (l®rj) + |*rj I^Jp^ (l^rj) 

ei €JV e ® 

+ l^rj • 


By using (29), we obtain 

v Vp k xp k < — ^xp k xp k = -VT k Vr k 

which implies that the xj - k -subsystem is ISS. 

Step 2: We define x Si = {xp k '■ e-k £ Qp and ek £ £' c } as the state of the strongly connected 
component. By taking (29) as the input, each subsystem xp k is ISS and admits an ISS-Lyapunov 
function Vp k . From Lemma 2, for the set of all the simple loops 67 -, if the cyclic-small-gain 
condition 


A tMZ) < Id 

is satisfied, then the composed subsystem x Si is ISS. By taking the strongly connected subsystems 
x Si as nodes, the upside subsystem Hp is acyclic as we previously mentioned in Lemma 8 . From 
[24], we note that the ISS properties are retained if the underlying digraph is acyclic. Therefore, the 
upside subsystem Hp is ISS as well. 

Additionally, we can verify that Vp defined in (27) is an ISS-Lyapunov function. Also we can 
calculate the the interconnection gain -/H from He to Hp. To begin with, according to 

Vn - { 7 S iVc) } 


we have 


Vp>(N- lfrg. (V c ),k = 1,2,--- ,N-l 

N-l 

since Vp = £ v n ■ 

k =1 

By choosing crk = cr > 0, we obtain 

N—l N-l 

\7Vpxp = X v v T k xp k < - X y x T k x T k = -crVp 

k=1 k=1 

which implies Vj- k is an ISS-Lyapunov function. Then we can simply choose 7 ~ c as 

l! c T = (N-lht Tk . (32) 

Step 3: Since we have xc = T T xp from (11), then 

v c = \xlx c = ^xpTT T x r 
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where TT 1 is symmetric positive semidefinite. Suppose the eigenvalues of TT t can be ordered and 
denoted as 0 < Ai < A 2 < • • • < \n-i- Clearly, one can obtain that Vc > Ai Vj-. 

Assume that P is an orthogonal transformation matrix and let xj- = Py, then we can translate Vc 
into a standard quadratic form as follows: 

Vc = + A2J/2 + • • • + A n-iVn-i)■ 

By taking the derivation of Vc, one can obtain 


N-l N -1 

W c £e = 22 Xkvlilk = 22 XkXrPkP^xr 

fc=l k- 1 

N-l N-l N-l 

= ^ A kWrxr < —cr XkYr < — =-A kVc 

k= 1 k =1 Xn ~ x k=i 

which implies that Vc is an ISS-Lyapunov function. Then we can choose the interconnection gain 
as 


Tf c r = Ai. (33) 

For this two interacting subsystems Hj- and He, if the small gain condition -V c o 7 /J < Id is hold, 
then the whole system is ISS. To satisfy the small gain condition, by combining (32) and (33), we 
can choose 


< Ai (N - 1) ‘ 


(34) 


From Theorem 1 , it is clear that the objective robust consensus can be guaranteed while (34) is hold. 

□ 


Remark 4 

The ISS-Lyapunov function for the composite system x (t) can be obtained by using the approach 
mentioned in [111. Besides, the discussion about the explicit cyclic-small-gain conditions required 
in (3) can be found in our previous study [26]. In particular, we can check that the cyclic-small-gain 
conditions can be guaranteed by simply choosing the nonlinear gains 7 - Ci < Id, with e/ £ N®. 

e k 


5. SIMULATIONS 

Numerical simulations are performed to illustrate the obtained theoretical results. For this set of 
simulations, we consider a six-agent system with both strongly connected graph and quasi-strongly 
connected graph. The dynamics of the 7-th agent is assumed to be 

Xi (t) = / (t, Xi ) + Wi (t) + Ui (t ), 7 = 1,2,--- ,6 

where 27 (t), m (t) £ R 3 , with the inherent nonlinear dynamics / : R x R 3 —> R 3 described by 
Chua’s circuit 


/ (t, Xi) = (C {-Xu +Xi2-l {xn)) , Xu - Xi2 + X i 3 ,-XXi2) T (35) 

where l {xu) = bxu + 0.5 (a — b) (\xa + 1| — \xu — 1|). While choosing ( = 10, x = 18, a = 
—4/3 and b = —3/4, system (35) is chaotic with the Lipschitz constant r\ = 4.3871 ([30]) as shown 
in Figure 8 . Assume the state of each agent is corrupted by white noise 117 (t) £ R 3 with the noise 
power £ = [0.25 0.25 0.25] T . 
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Figure 8 . Chaos with two attractors. 


5.1. Case 1: Strongly Connected 

The digraph is strongly connected as shown in Figure 5(a). From (22), multi-agent system can 
be translated into the edge Laplacian dynamics associated with the edge-interconnection digraph 
shown in Figure 5(b). The incidence matrix E and the edge adjacency matrix A e are 


E = 


Ae = 


/- 1 

1 

0 

0 

0 

0 

0 

°\ 

1 

0 

-1 

1 

-1 

0 

0 

0 

0 

-1 

0 

0 

0 

1 

-1 

1 

0 

0 

1 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

0 

0 

V 0 

0 

0 

0 

0 

0 

1 

-1/ 

/° 

-1 

0 

1 

0 

0 

0 

°\ 

0 

0 

0 

0 

0 

-1 

0 

-1 

-1 

0 

0 

-1 

0 

0 

0 

0 

1 

0 

-1 

0 

0 

0 

0 

0 

-1 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

1 

0 

0 

0 

0 

0 

-1 

0 

-1 

V 0 

0 

0 

0 

0 

1 

-1 

0 / 

be calculated though Le 

= I 

+ Ae. 




By simply choosing 7 ^' < Id with fc, Z = 1,2,--- , 8 , e; S IV®, the cyclic-small-gain theorem 
condition is satisfied. By taking (s) = 0.9487s and a = a = 1, then we obtain 


X e , _1 

Pi" =QL o 


o a (s) = 1.0541s. 


After choosing a ek = 1, the input for the edge-interconnection system (20) is proposed as 


ttei — 1 


Ue , = 




V \X e 


2.1082 \x e 


2.1082 lie 


+0.5 x ei ,u e2 =-- 
Vzt^J +0.5 Xe 3 ,U ei = -- 


V \Xe 2 \ 

v\x ei \ 


2.1082 lx e 


2.1082 \x e 


+ 0.5 Xe 5 ,U ee = — 


T]\x e6 \ + 2.1082 |m e6 1 + 


^2f) 

vk) 

vk) 


U e7 = —; 


0.5i e 


0.5i e 


0.5i P 


V \x e5 1 + 2.1082 \x et \ + v / 2^) 

(j] \x e7 \ + 2.1082 \x er \ + V2^j + 0.5i e7 , it es = (v |^e 8 | + 2.1082 |x es | + + 0.5i e 


Finally, by using (24), the consensus protocol ( 14) can be obtained. 
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The simulation results are shown in Figure 9. The edge states x ek converge to the neighbors of 
the origin by applying the consensus protocol shown in Figure 9(a). From Figure 9(b), 9(c) and 
9(d), one can see that the robust consensus is indeed achieved. Therefore, the proposed consensus 
protocol can effectively address the challenges resulting from the inherently nonlinear dynamics 
and unknown but bounded disturbances. 




(a) The edge agreement is achieved under control law (b) The evolutions of xn. 

(23) and (24). 




(c) The evolutions of (d) The evolutions of xa. 

Figure 9. The simulation results under strongly connected digraph. 


5.2. Case 2: Quasi-strongly Connected 

In this case, Figure 6(a) depicts a quasi-strongly connected digraph, while Figure 6(b) illustrates 
its corresponding edge-interconnection digraph. The yellow nodes in Figure 6(b) correspond to Qy. 
According to the partition, the incidence matrix of the spanning tree Ef, the incidence matrix of the 
cospanning tree Eq and the edge-adjacency matrix A e are as follows: 


Ej- 


ei 

e-2 

e 3 

e 4 

67 


e5 

e6 

( 1 

1 

1 

0 

0 ^ 


( 0 

° \ 

-1 

0 

0 

1 

0 


-1 

0 

0 

0 

-1 

0 

1 

Ec = 

0 

1 

0 

0 

0 

-1 

0 

0 

0 

0 

-1 

0 

0 

0 


1 

-1 


. : -1/ V 0 0 / 
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ei 

e-2 

£3 

e 4 

e7 

es 

ee 

ei 

(° 

1 

1 

-1 

0 

0 

° \ 


1 

0 

1 

0 

0 

-1 

0 

63 

1 

1 

0 

0 

-1 

0 

-1 

e 4 

0 

0 

0 

0 

0 

0 

0 

e-i 

0 

0 

0 

0 

0 

0 

1 

£5 

0 

0 

0 

-1 

0 

0 

0 

ee 

V 0 

0 

0 

0 

1 

-1 

0 / 


By simply choosing 7 ^ (s) = 0.9487 s,k,l = 1,2,3,4,7, e; € iV®, the cyclic-small-gain condi¬ 
tion is satisfied. From (12), we could have 


T = 


' 1 

-1 

0 

0 

. 0 


0 ■ 
1 

-1 

0 

0 . 


The smallest non-zero eigenvalue of TT t is Ai = 1; therefore, from (30), we can choose 77 ^ = 
0.175s and resulting 

p7t, = “ _1 0 ( 7 ^' ) 0 a (s) = 1.0541s 
pZ k = a ” 1 0 ( 7 t c Tk ) 0 a (s) = 5.7143s. 

After taking a ek = 1, the control input for each edge-interconnection system (20) can be determined 
as 


u ei 

u e 3 

U ei 

Ue 7 


XT l 
l®Til 

(v\XTi\ 

+ 3.1623 |x Tl | 

+ 

+ 0 . 5 x 77 

*r 2 

\xt 2 \ 

(v \xt 2 \ 

+ 7.8225 \xt 2 1 

+ V2£) 

+ 0 . 5 x 77 

xt 3 

\xt 3 \ 

( v\xt 3 \ 

+ 8.8766 \xt 3 1 

+ V2Z) 

+ 0 . 5 x 77 

XTi 

\xTi\ 

(v\xtA 

+ + 0 . 5 x 77 


X Tl 

\xr 7 \ 

(v\xt 7 1 

+ 5.7143 \x T7 \ 

+ V2£) 

+ 0 . 5 x 77 . 


Finally, based on (14) and (31), the implementable consensus protocol can be obtained. 

Figure 10(a) shows that the edge states reach agreement by using the consensus protocol. 
Simultaneously, multi-agent system achieves robust consensus shown in Figure 10(b), 10(c) and 
10(d). Clearly, the proposed consensus protocol can effectively restrain the influences resulting from 
the inherently nonlinear dynamics and the unknown but bounded disturbances. 


6 . CONCLUSIONS 

The edge Laplacian of digraph and its related concepts were originally proposed in this paper. 
Based on these graph-theoretic tools, we developed a new systematic framework to study multi¬ 
agent system in the context of the edge agreement. To show how the edge Laplacian sheds a new 
light on the leaderless consensus problem, the technical challenges caused by the unknown but 
bounded disturbances and the inherently nonlinear dynamics were considered; and the classical 
ISS nonlinear control methods together with the recently developed cyclic-small-gain theorem were 
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(a) The edge agreement is achieved under control (b) The evolutions of am. 

law (29) and (31). 




(c) The evolutions of am . (d) The evolutions of am . 

Figure 10. The simulation results under a quasi-strongly connected digraph. 


successfully implemented to drive multi-agent system to reach robust consensus. Furthermore, the 
edge-interconnection graph, which plays an important role in the analysis and synthesis of multi¬ 
agent networks, was proposed and its intricate relationship with the original graph was discussed. 
For the quasi-strongly connected case, we also pointed out a reduced order modeling for the edge 
agreement in terms of the spanning tree subgraph. Based on this observation, by guaranteeing the 
ISS properties of each subsystem and assigning the appropriate gains for both of the interconnected 
subsystems to satisfy the small gain condition, the closed-loop multi-agent system could reach 
robust consensus. Under the edge agreement framework, we believe nonlinear multi-agent systems 
with more complex factors, such as switching topologies and time-delays, will be well settled. 
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